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Abstract 


A class of Bernstein-like basis functions, equipped with a shape param- 
eter, is presented. Employing the introduced basis functions, the corre- 
sponding curve and surface in rectangular patches are defined based on 
some control points. It is verified that the new curve and surface have 
most properties of the classical Bézier curves and surfaces. The shape 
parameter helps to adjust the shape of the curve and surface while the 
control points are fixed. We prove that the proposed Bézier-like curves 
can preserve monotonicity and that Bézier-like surfaces can preserve axial 
monotonicity. Moreover, the presented curves and surfaces preserve bound 
constraints implied by the original data. 
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1 Introduction 


In computer-aided geometric design (CAGD) and computer graphics, curves 
are often constructed by the following relation: 
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CH => AOR t € [a,b] CR, neEN, (1) 


where the real functions {F;} are known as blending functions, {P;} C R* 
are control points (s > 2), and the polyline constructed by control points is 
called the control polygon. One of the most important blending functions 
are the Bernstein polynomials, which in turn lead to well-known classical 
Bézier curves. They have a simple structure and are widely used in many 
engineering and technology fields. 

In a classical Bézier curve, once the control points are fixed, the shape 
of the curve cannot be changed. This is a drawback that reduces the shape 
adjustment and therefore limits their applicability. In order to overcome this 
deficiency, many researchers have tried to construct basis functions equipped 
with shape parameter(s) to create new curves whose structures and properties 
are similar to the Bézier curves, besides they also have the shape adjustment 
property [6, 13, 16, 17, 19, 21, 22, 23, 30] 

In curves with shape parameter(s), by changing the parameter(s), the 
curve either approaches to the control polygon [17, 18, 19], or moves away 
from it [20, 31]. Another difference between parameter based curves origi- 
nates from their construction, some basis functions are defined by a recur- 
sive formula [23, 22, 5, 20, 32], while others are based on a general formula 
[17, 19, 16, 18, 24, 8]. 

Here, we propose a new family of parameter based blending functions 
by using square roots of polynomials. Thus we call them sq-basis functions. 
The new blending functions are generated from three initial basis functions 
employing a recursive relation. This is a procedure that could be traced back 
in the literature [5, 22, 23, 32]. The corresponding Bézier-like curves, which in 
turn are referred to sq-curves, are in common with the classical Bézier curves 
in many features such as geometric invariance, endpoint interpolation, formal 
symmetry, and convex hull property. The sq-basis functions are equipped 
with one shape parameter, which can adjust the shape of the curve, and they 
would either move closer to or away from the control polygon. In the sq- 
curves with three control points, by altering the shape parameter, one would 
see a family of curves traveling from the control polygon to the straight line 
joining the first and last control point. In comparison to previously defined 
blending functions with the same structure, sq-curves show a higher rate of 
change by employing just one shape parameter, and this could be considered 
as an advantage. The construction could be extended to surfaces [11, 14], 
which share the most common features of Bézier surfaces. 

Once we have a suitable model for data approximation, one can apply it 
in a shape-preserving approximation. Shape-preserving approximation is de- 
fined to be a method of constructing a curve (surface), which also preserves 
the shape implied by the data. It is known as an essential curve/surface 
design technique in CAD/CAM and geometric design. Convexity, mono- 
tonicity, positivity, and boundedness are the most important shape features, 
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which have extensively been studied in literature [1, 2, 3, 25, 27, 28, 29]. The 
proposed Bézier-like curves and surfaces prove to be a promising approxima- 
tion tool that preserves the monotonicity and the boundedness implied by 
the data set. 

The present paper is outlined as follows. In Section 2, the new blending 
functions are defined and their properties are studied. The corresponding 
curves and surfaces are defined in Section 3. Section 4 studies the shape- 
preserving properties of the new curves and surfaces. 


2 New blending functions 


The classical Bernstein functions have a number of important features that 
makes them a suitable basis for constructing control point based curves. Non- 
negativity, partition of unity, symmetry in some sense, and their end point 
values are among the most cited properties of classical Bernstein functions. 

Here we address a general representation for a class of blending functions 
that have similar properties to Bernstein polynomials. We present a special 
case that uses the square root of quadratic (sq) polynomials. 


Definition 1. Based on a real valued shape parameter v, the sq-starting 
basis functions are defined as follows: 


bolt) = 5 -t + 0), 
bn. (t) =1-29(t), (2) 
bolt) =t- 5 + (8), 


We generate the sq-basis functions of order n (n > 3) with the recursion 
relation of the classical Bernstein basis functions [9] as follows: 


bri (t) = (1 a t)bn—1,i(t) + tbn—1,i-1(t), te [0, 1], (3) 


where i = 0,1,2,...,n. For k = —lork>n, we set bp,(t) = 0. 
The sq-basis functions of order n = 3 can be expressed as 


bolt) =P -St+ 5+ (1-H eld, 
b3.1(t) = —¢? it + 1+ (—2+ 3t) y(t), (4) 


IJNAO, Vol. 12, No. 2, pp 449-466 


452 Nouri and Saeidian 


n=2 = 
1 1 + r + : : 1 bcd 
a 0.01 
09 WN 09 0.1 
— y=0.3 
08 ~S v=0.6 
S oe —— 0.8 
“ —— ml 
07 ~ 07 
06 06 
y 
05 >< ost - me 
04 S 04 
03 03 
02 a: 02 - eS. 
ee SS 
o4 - one o4 EZ SS 
é Zo . 
le ie Hes —— = 
9 L= a =— —=> 0 i = 
0 o1 02 03 04 05 06 O07 08 09 1 0 01 02 03 04 O05 06 O7 08 09 1 
n=4 n=5 
1 + 1 
w=0,01 v=0.01 
ost v=0.1 09 v0.4 
w=0.3 v=0.3 
ie v=0.6 08 10.6 
98 v=0.8 0.8 
sia wt 
07 07 
06 06 
05 05 
04 _— L 04 - — 
03 SS "ear SOY, 03 —s 
< OS EY, = 
02 tty, 0.2 
04 = orf Zz 
= = = ff = = 


Figure 1: The sq-basis functions with different values of v. 


b3,o(t) = -t* + <t- 5 + (1— 3t) y(t), 


1 
b3 3(t) = — at + ty(t). 


Figure 1 shows the sq-basis functions generated with n = 2,3,4,5 and differ- 
ent values of v. 


Theorem 1. The sq-basis functions generated by (2) and (3) have the fol- 
lowing properties: 


(a) Nonnegativity: b,;(¢) > 0 for i =0,1,2,...,n. 
(b) Partition of unity (Normalization): )7}" 9 bn,i(t) = 1. 
(c) Symmetry: bpi(t) = bnn—i(1 — t) for i = 0,1,2,...,n. 


(d) End-point values: 
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—n+v, 7=0, n-V tn, 
n—-2v, t=1, 2av—n, t=n-—1, 

COS ee ee Se) 
0, 0.W., 0, 0.W. 


Proof. Mathematical induction is used to prove this theorem. 


(a) Nonnegativity is obvious from (2) and (3). 


(b) 


From (2), we obtain per b2,i(t) = 1. Now suppose that the equality is 
true for m. Then according to recurrence relation(3) and the inductive 
hypothesis, we obtain 


m+1 m+1 m+1 
S> bm+1i(t) = (1—t) 35 ome) +t 55 Om sit) =1-t+t=1. 
1=0 1=0 1=0 


The sq-starting basis functions are symmetrical. Assume that the sq- 
basis functions of order m are symmetrical. Then from this inductive 
hypothesis and recurrence relation (3), one has 


bm+i1,(1 —#) = (1 — (1 —t))bm (1 — t) + (1 — t)bm s_1(1 — £) 
= tbin,m—i(t) + (1 a t) bm, m—i+1(t) = bm+1,m—i41(t). 


From (2), we have 


where 


By a simple deduction from (2) and (7), we conclude that the results 
in (5) and (6) hold for n = 2. 


Suppose that the properties at the endpoints hold for the sq-basis func- 
tions of order m. Then 

1 i=0 

114(0) = Pm (0) i ag 


1 t=m, (8) 


bm+1,i(1) = bmi (1) = 0, i#m 


which can be obtained by the inductive hypothesis and (3). 
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From (3), we have 


Bing1,i(0) = —bm,i(0) + Bin,i(0) + bm,i—-1(0), (9) 
Bin g1,4(1) = —Omng (1) + be (1) + Ong )- (10) 


From (8) and the inductive hypothesis and by using relation (9), the 
following results are obtained: 
(i) For ¢ = 0, bin41,0(0) = —bm,o(0) + On,0(0) = —1- (-m+y) = 
= (m + 1) Ty; 
(ii) For i = 1, bi,411(0) = —bm,1(0) +8), 1 (0) +Om,0(0) = (m + 1)—2y, 
(iii) For 4 = 2, bin419(0) = —bm,2(0) +0), 2(0)+0m,1(0) = O+v+0 = », 
(iv) For i £0,1,2, 6,41 ;(0) =0. 


Similarly, from (10), the following results are deduced: 


(i) Fori=m-+1, Oia!) =(m+1)-y, 
(ii) For 7 =m, bf (1) = —(m+1) + 2p, 


m+1,m 


) 
(iii) For i= m—1, Unsim—1(1) = —-», 
) 


Mm 


(iv) Forti Am—1,m,m+1, bai i(1) = 0. 


3 sq-curves and sq-surfaces 


The sq-basis functions can be considered as blending functions to construct 
Bézier-like curves. Here we present the sq-curves and also extend the idea to 
surfaces. 


Definition 2. Given control points {P;}i_, CR’, for0<t<1,0<v<1, 


rn(t) = > bnilt)Pi, (11) 


is called a sq-curve, where {b,,;(t)}\.5 


defined in (2) and (3). 


(n > 3) are the sq-basis functions 


From Theorem 1, the sq-curve introduced in (11) has the following prop- 
erties: 


(a) Geometric property at the endpoints: From the properties at the end- 
points of the sq-basis functions, we obtain 
rm(0) = Po, rn(1) = Ph, 
ri (0) = (v) (P2 2P, t Po) t n(Pi Po), 
ri (1) = (-v) (Pa — 2Pr—1t+ Pa—2) +7 (Pn - Pa—1)- 
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Figure 2: sq-curve with different values of v. 


(b) Symmetry: The symmetry property of the sq-basis functions leads us 


to 
Mn(t; Po, Pi, Pa sich) = So bn i(t)Pi = Peon el = t)P; 
1=0 1=0 
= pa e! — ee — rn(l = te Pais: . ., P). 
j=0 


(c) Geometric invariance: As r,,(t) is an affine combination of the control 
points, so the shape of the sq-curve is independent of the choice of 
coordinate system. 


(d) Convex hull property: Because of the nonnegativity and normalization 
properties of the sq-basis functions, one concludes that the sq-curve 
must be inside the convex hull of its control polygon. 


Figure 2 shows two specific curves generated by sq-curves, in which the 
control points are fixed and values of v are set to v = 0.01, 0.2, 0.4, 0.6, 0.8, 1, 
respectively, from outside to inside. One can see that the sq-curves are to 
approach the control polygon as the shape parameter decreases. 


3.1 Comparison with existing basis functions 


The construction of a suitable set of blending function has attracted a great 
deal of interest among scientists. The main goal of most papers in this filed 
is to construct parameter based blending functions, which by altering this 
parameter(s), one can make most changes in the corresponding parametric 
curve (1). Some works focus on a specific application, and others studied 
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the structure and construction strategies. Here, we aim to emphasize on 
advantages of the sq-basis functions. As the sq-basis functions are generated 
recursively, we compare them with those who have a recursion formula. 

Figure 3 represents a visual comparison. Figure 3(a) shows the sq-curves; 
Figure 3(b) represents the curves constructed in [32]; and the curves proposed 
in [23] are depicted in Figures 3(c) and 3(d), for different values of shape 
parameters. Compared to [32], the rate of change in sq-curves is greater than 
changing the corresponding shape parameter. The curves generated in [23] 
show a good rate of change for different values of shape parameters, but this 
is achieved by employing two separate parameters. In the sq-curves, this 
much of flexibility is gained by altering just one shape parameter. 

In [20], no curve can be constructed for three control points because the 
starting basis functions have four members. Authors in [5, 6, 22, 23] employed 
three starting basis functions, but several shape parameters were used to 
make further changes to the curves. However, in the case of sq-curves, the 
same changes are gained with just one parameter. In some works [31, 15, 4], 
only four [15, 4] or five [31] bases were presented, and curves with fewer 
control points cannot be represented by these basis functions. For more 
control points, the continuity conditions are checked. 

Generally, among the Bézier-like curves, which aim to present curves that 
move from the convex combination of the first and last control points to the 
control polygon by altering shape parameter, the sq-curves serve as a good 
choice and compete well with the other bases. 


3.2 sq-surface 


The corresponding sq-surfaces can be constructed by a tensor product ap- 
proach. Given control points {P;};_) C R®, the surface represented by 


S(8,t)=> >>) Pizbmi (a) ong (), O<8,t<1, (12) 


i=0 j=0 


is called a sq-surface, in which b,,,;(s) and b,;(t) are the sq-basis functions. 
Figure 4 shows that different surfaces can be obtained by altering the 
shape parameter when the control net is fixed. 
From the properties of the sq-basis functions, it is easy to verify the 
following properties of the corresponding sq-surfaces: 


(a) Interpolation at the corner points: S (s,t) passes through corner points 
Poo, Pm,o, Pon; and Pm». In fact, we have $ (0,0) = Poo, S$ (1,0) = 
Pm,o, 5 (0,1) = Pom, and S$ (1,1) = Pinjn- 


(b) Nonnegativity: It is obvious that bm; (s) bn,; (t) is nonnegative for all 
m,n,i,j and s,t € [0,1]. 


IJNAO, Vol. 12, No. 2, pp 449-466 


On a class of Bézier-like model for shape-preserving ... 457 


(c) v=—-1, p= —1,-0.6,—-0.2,0.2,0.6,1 (d)v=1, »=—1,-0.6,—0.2,0.2,0.6,1 


Figure 3: Comparison of sq-curves with curves proposed in [32] and [23] 


Figure 4: Shape of sq-surfaces with different shape parameters. 
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(c) Partition of unity: 


>: 3 bmi (s) bn j (t) = 3 bmi (s) 3 bn, j (t) = 1, 0 < St < 1. 
i=0 7=0 1=0 j=0 


Proof of this is due to the partition of unity of the sq-basis functions. 


(d) Convex hull property: The surface lies in the convex hull of its control 
points, since S (s,t) is the linear combination of its control points with 
positive coefficients that sums to 1 (partition of unity). 


(e) Affine invariance: This means that the surface defined by control points 
that are an affine transformation of all control points is the same surface 
that is obtained by applying the same transformation to the surface’s 
equation. As S(s,t) is an affine combination of the control points, so 
the shape of the sq-surfaces is independent of the choice of coordinate 
system. 


(f) Boundaries of surface: The four surface boundaries are exactly the four 
sq-curves defined by the boundary edges of the surface as follows: 


j=0 g=0 
S (8,1) = 5° Pinbmi(s), $ (8,0) = 5° Piobm, (8) 
i=0 i=0 


4 Shape-preserving approximation by sq-curves and 
sq-surfaces 


The sq-curves preserve monotonicity and the sq-surfaces are axially mono- 
tonicity preserving. The newly defined models preserve boundedness in both 
two-dimensional and three-dimensional data. Here we state these facts in a 
detailed discussion. 


4.1 Monotonicity preservation of sq-curves 


Here, we present a proof of the monotonicity preservation of the sq-basis 
functions and sq-curves. 


Definition 3. A system of functions (go,..., gn) is monotonicity preserving 
if for any Ag < Ay <-+: < An in R, the function )>"_9 Aig; is increasing. 
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The following result, which characterizes monotonicity preserving sys- 
tems, appears in [7, Proposition 2.3]. 
Proposition 1. Let (go,...,gn) be a system of functions defined on an 
interval [a,b]. Let vj := )(j_,9; for i € {0,1,...,n}. Then (go,-.-,9n) 
is monotonicity preserving if and only if vo is a constant function and the 
functions v; are increasing for 7 = 1,...,n. 

We obtain the following result using the previous proposition. 


Proposition 2. The sq-basis functions (bn.o,bn1,.--,0n.n) defined by (2) 
and (3) are monotonicity preserving. 


Proof. Let us consider sq-basis functions (2) with n = 2. Using Proposition 
1, (62,0, b2,1, 62,2) is monotonicity preserving if vp is a constant function and 
the functions v; and v2 are increasing. Moreover, 


2 
Ug = S/ be,5(t) = 1 
j=0 


and 
11 = Do baglt) = baa) + boa) = 5 +t- lt) + vil) = 1-90. 


The function y(t) is a convex function in the interval [0, 1], so the decreasing 
function —y’(t) takes its minimum at ¢ = 1. Thus we have v{(1) = v, so 
v(t) > 0 for allt € [0,1]. Therefore the function v(t) is increasing. We 
have 


Again, employing the convexity of y(t), one observes that the increasing 
function y’(t) takes its minimum at t = 0. Since v5(0) = v results in v(t) > 
0, for all t € [0,1], so we can conclude that v2(t) is an increasing function. 


Suppose that the sq-basis functions (b,—1,0, bn—1,1,- ++; 0n—1,;n—1) is mono- 
tonicity preserving. Then using Proposition 1, (6n,0,bn1,---,0njn) is mono- 
tonicity preserving if vg is a constant function and the functions v; are in- 
creasing fori =1,...,n. 

We have 


vo = So bn (t) = 1 
j=0 
and 


1 = Yo bnglt) > Vi(t) = S705 (0 a 
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By combining formulas (3) and (13), we have 
=o ig(t +09 oh 1,9 (t)+ So bn—15—1( 1 Oth 1jj- i(t) 
j=i 


n-1 n-1 
= bn-in+|—- S- bn—1,5(t) + S- bn—1,5(t) | + bn—1,1-1 (4) 
ji ji 


n 


+(1-t) da Lj OF Ed tr-rgall 


= bp_14-1(8) +09, ia Sue ie ag EY (14) 


j=i-1 


Since (bn—1,0, 0n—1,1,---,0n—1,n—1) is monotonicity preserving, we deduce 


yh. ig) 2.0, yb () \ 0. (15) 


j=i-1 


Now, from (14) and (15), for ¢ € [0, 1], it is obvious that vj(t) > 0. 


The sq-curve is a curve in R? plane, so we need the following statement 
to assure the monotonicity preservation of sq-curves. The proof is straight- 
forward. 


Proposition 3. Let (go,..., gn) be a system of functions that is monotonic- 
ity preserving. If {P;};) C R? is the monotone data, then the function 
eo Pig: is increasing. 


4.2 Monotonicity preservation of sq-surfaces 


The sq-surfaces (12), defined on rectangular patches, are axially monotonicity- 
preserving. To verify this assertion, we state our reasoning based on some 
known results from the literature. 

Given two strictly increasing sequences of abscissas aw = (a, Q1,..-,Q@m) 

i p=Osis0570- 

and 6 = (89, 81,...,8n) and control points {(2 8; i} me , we define the 
control net p : [ao, @m] X [80, Bn] > R to be the unique function that satisfies 
the interpolation conditions p(a;,6;) = cj, for all i = 0,1,...,m, and j = 
0,1,...,m, and is bilinear on each rectangle R;; = [a;, ai41] (6), Bj41]- 

A bivariate function g is increasing in a direction d = (d1,d2) € R?, if 
g(x@t+Adi,y + Adz) > g(x,y), > 0. In particular, the control net p can be 


IJNAO, Vol. 12, No. 2, pp 449-466 


On a class of Bézier-like model for shape-preserving ... 461 


increasing in a direction d. This case has been characterized by Floater and 
Pefia [10] in the following way. 


Lemma 1. The control net p is increasing in the direction d = (di, dz) € R? 
if and only if for 7 =0,1,...,m—1 and j =0,1,...,n—1, 
dy Ai ci, 541 + dg A2ci+k,j > 0, k,l € {0, 1} ; 


where Aycj; := (ci+1,5 — Ci,j) / (iti — a4) and Ageys == (c,541 — i,5) / (A541 — B5)- 
O<j<n aa 
O<icm We have Aiciy : Cizi,j — Cy for i = 
.,m—1 and Jj = 0,1,...,n — 1 and Noci; : Cij4+l — Cig fori = 
.,m—landj=0,1,...,n—1. 


Given a sequence (c;;) 
0, ty. 
0,1,.. 


Remark 1. As a consequence of Lemma 1, we have that the control net p is 
increasing in the X-axis direction d = (1,0) if and only if for i =0,1,...,m-—1 
and j = 0,1,...,n—1, Aicj; > 0. Analogously, the control net p is increasing 
in the Y-axis direction d = (0,1) if and only if for i = 0,1,...,m—1 and 
j3=0,1,...,n-1, Noci; > 0. 


Theorem 2. If the control net p is increasing in the X-axis direction d = 
(1,0) or in the Y-axis direction d = (0,1), then so is sq-surface S (s,t). 


Proof. Suppose that the control net p is increasing in the X-axis direction 
d = (1,0). Then for 7 = 0,1,...,m—1 and j = 0,1,...,n —1, we have 


Aci; > 0. We show that, in this case, 5 (s,t) is increasing in the direction 
d = (1,0). Therefore, 


$= > 0 | 6 | ome @) ong, OX 94<1. 


The surface S(s,t) in the direction d = (1,0) for specified y = y is a curve 
that is obtained from the collision of the plane y = 7 and the surface S (s,t), 
se) 


S> >- Bjbmné (8) Ong (t) = 7 > J bmi (8) Biba, () = 
1=0 7=0 i=0 7=0 


=> >> Bjbng ®) =7- 
j=0 


Therefore ¢ is fixed and the curve on plane y = 7 is as follows: 
> oes (t) » (**) Daraci 2) , O<s<l. (16) 
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Now from Proposition 3, one observes that >)" 9 (**) bmi (8) is increas- 
ij 


ing for any j, which lead us to the conclusion that the curve (16) would be 
increasing. Because the value of 7 was arbitrary, the surface S$ (s,t) is in- 
creasing in the direction d = (1,0). One can prove the desired result in the 
direction d = (0,1), in a similar way. 


Example 1. The set of data presented in Table 1 is generated from the 
function 
F(2,y) =In(a? +42). (17) 


Ti j=O,...,6 
We have {( 5 )} as a data set in a three-dimensional space. 
F(ais93)/ J j=0,...,6 


These provide us 49 control points. By these control points, we have a 
control net p, which is increasing in the X-axis direction d = (1,0) and in 
the Y-axis direction d = (0,1). Figure 5 is the visual model of monotone 
surfaces in the X-axis direction and in the Y-axis direction that are obtained 
by control points and (12) with different values of free parameter v. The 
error rate of monotone surfaces has been reported. 


Table 1: Data generated from the function F(x, y) = In(x? + y?) 
y x 
20 40 60 80 100 200 300 

20 6.6846 7.6009 8.2940 8.8247 9.2496 10.6066 11.4120 
40 7.6009 8.0709 8.5564 8.9872 9.3588 10.6359 11.4252 
60 8.2940 =8.5564— 8.8818 = 9.2103 9.5178 ~=—:10.6828 =—:11.4468 
80 8.8247) = 8.9872 9.2103 9.4572 9.7050 =—-:10.7451 11.4763 
100 «9.2496 9.3588 9.5178 9.7050 9.9035 10.8198 11.5129 
200 10.6066 10.6359 10.6828 10.7451 10.8198 11.2898 11.7753 
300 11.4120 11.4252 11.4468 11.4763 11.5129 11.7753 12.1007 


4.3 Bound preservation of the sq-model 


Suppose that we have a discrete data that is generated from a bounded 
phenomenon, either in a two-dimensional space or in dimension three. When 
we attempt to visualize the data using approximation methods, it is natural 
to expect the approximation to adhere to the boundedness [26, 25]. The sq- 
curves and sq-surfaces fulfill the convex hull property, which is the powerful 
feature that guarantees the bound-preserving approximation. Actually this 
is a feature that is in common for all Bézier-like model with convex hull 
property. When we approximate a data with a Bézier-like model, the curve 
or surface lies in the convex hull of the control points. If we set the data to 
be the control points, then the curve (surface) would be entirely inside the 
control polygon (control net), which guarantees the boundedness. 
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Figure 5: Axial monotonicity-preserving approximation with different shape 
parameters and the error functions 
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5 Conclusion 


We presented Bézier-like parametric curves and surfaces through introducing 
a new class of blending functions with one shape parameter. These curves and 
surfaces, in addition to covering many properties of classical Bézier models, 
can be adjusted in a shape by altering the shape parameter when the control 
points are fixed. It is been verified that the newly defined sq-models can 
be successfully used for monotonicity-preserving approximation. The bound 
preservation property of the sq-models has also been reviewed. 

Note that the introduction of shape parameter not only brings an ad- 
vantage to the shape modification of the Bézier-like curves but also provides 
an optimized parameter for the shape optimization design of the Bézier-like 
ones. Hence, the research on how to utilize the SDABWO [12] to solve the 
model of curve shape optimization, which takes the shape parameter as the 
optimization variables, could be addressed in future works. 
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